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Abstract 

Non-Hermitian VT— symmetric Hamiltonians H = —d 2 /dx 2 + V(x) with 
x G K are reinterpreted as describing the most elementary phenomenological 
quantum graph, i.e., a system living on the two half-line edges connected 
at a single matching-point vertex in the origin. A q— pointed star graph 
generalization of these q = 2 models is then proposed and studied. For a 
special toy-model point interaction yielding the exactly solvable model at 
q = 2, the bound-state energies are finally identified with the roots of a 
remarkably compact trigonometric function at any q — 2, 3, 
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1 Introduction 



The heuristic use of the concept of VT symmetry, i.e., of the parity-times- 
time-reversal symmetry of quantum Hamiltonians H and/or of the toy- model 
wave functions ip(x) (with, say, x G M) proved unexpectedly productive in 
phenomenologically oriented quantum field theory [TJ [2] or in the context of 
relativistic quantum mechanics [3]II], in the supersymmetric model-building 
[HIE] or, recently, in experimental classical optics [TJ. 

One of the simplest illustrative examples of a VT symmetric Hamiltonian 
has been proposed in Ref. [8]. In the model the quantum motion remained 
free inside a finite interval of coordinates, 

-^!(x) = E*(x), xe(-L,L). (1) 

The only dynamical information has been carried by the very specific point 
interaction induced by the external Robin-type boundary conditions contain- 
ing the single real coupling constant a, 

*'(±L) + ia *(±L) = , a > . (2) 

This dynamical input yielded the phenomenology-oriented real bound-state 
spectrum 

E = a 2 , E n = (-) , n = l,2,... . (3) 

Its explicit form enabled us to restrict our attention, for the sake of sim- 
plicity, to the non-degenerate systems where 2La/ii ^ 1,2, . . .. The proba- 
bilistic quantum-mechanical interpretation of the closed-form wave functions 
of the model also appeared feasible. The explicit formulae yielding the uni- 
tary forms of the model (i.e., in the notation of review j9], all of the non- 
equivalent "standard" inner-product representations of the eligible physical 
Hilbert space of states V, {s) ) were found and described in a series of mathe- 
matically rigorous subsequent studies [TUt ITTj. 

The appealing and, in some sense, extreme simplicity of the latter ex- 
ceptional model inspired, naturally, a number of generalizations (cf., e.g., 
Refs. [HI [T2] ). Also our present paper will describe a new generalization of 
the model. 
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2 A quantum-graph reinterpretation of the 
VT— symmetric square-well models 



A formal core of our present considerations will lie in the reinterpretation of 
Eq. (JTJ where the interval of x E (-L, L) will be treated as a union of a 
pair of equal-length subintervals (or "edges") e + = (0, L) and e_ = (— L, 0) 
forming an elementary "graph" with the single "vertex" at x — 0. 

In such it is necessary to distinguish between the theoretical and 

purely phenomenological informal aspects of such a reinterpretation. In- 
deed, the latter, "realistic" aspect is very natural. Traditionally, it finds its 
widespread use in quantum chemistry where, typically, the valence electron 
of an organic molecule may be often treated as moving just strictly along the 
atomic-bond edges [13] . 

In the former, more abstract and less phenomenological setting the restric- 
tion of the motion to the edges of a suitable graph may enormously simplify 
the underlying Schrodinger equation [H]. Recently, this idea made the study 
of various quantum-graph models extremely popular. Pars pro toto the in- 
terested reader may be recommended to consult a comprehensive collection 
[15] of more than 700 pages of reviews and original research reports, with the 
scope ranging from certain entirely "unrealistic" scenarios (i.e., e.g., from 
the fractal and/or chaos-simulating graphs [16] ) down to certain very realis- 
tic models of observable photonic crystals and various other "leaky-graph" 
nanostructures encountered, typically, in condensed matter physics [17] . 

For the sake of definiteness, let us now assume that in our above most 
elementary graph G^ 2 \ both of the respective edges are oriented inwards, 
i.e., e± = e±(y±) with y + = L — x 6 (0, L) while y_ = L + x 6 (0,L). 
Without any real loss of generality, our attention will also remain restricted 
to the dynamics represented by the end-point point interaction as mediated 
by boundary conditions (|2J). 

In the new notation we have to replace Eq. (pQ) by the pair of differential 
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Schrodinger equations 

-■frMv±) = EMv±), v±e(P,L) (4) 
a y± 

complemented by the standard regular matching conditions in the origin, 

ML) = ML) , d y+ ML) + d y _ML) = . (5) 

The external Robin-type boundary conditions © must now read, mutatis 
mutandis, 

V±(0) =±iaV±(0). (6) 

The physics (i.e., the spectrum) remains unchanged but the mathematical 
meaning of the VT— symmetry of H (i.e., the representation of the antilinear 
operator uj = VT) gets modified. 

2.1 Operators of symmetries 

After the change of the language, the differential-operator Hamiltonian (orig- 
inally defined as acting, in general, in L 2 (M)) must be treated as acting in 
another, "friendly" [9.J Hilbert space of states U {F) = L 2 (M + )®L 2 (R+). 
Formally, Schrodinger Eq. (j4]) then acquires the two-by-two operator-matrix 
form 

H + — E 
H_ -E 
(i.e., (H — EI)\ifj) = in an abbreviated notation). 

The original linear operator of parity V (i.e., the reflection which changed 
the sign of the coordinate, V : G (—L,L)) will now play the 

slightly different role of a domain-intertwiner such that V : y± — > y T , i.e., 

<8> 

The parallel interpretation of the antilinear symmetry u> = VT of H will vary 
with the spectral properties of H [T5]. As long as we may define \ip') = u)\ijj), 
we may write u(H — EI)\ip) = (H — E*I)\ip') = 0. Thus, in the generic 
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non-degenerate case we must distinguish between the real-energy scenario 
E = E* (in which will be proportional to and the case of E ^ E* 
in which the two eigenvectors and of H remain linearly independent 
(or, in the language of wave functions, in which the VT symmetry becomes 
spontaneously broken [2]). 

2.2 Wave functions and energies 

In the case of the unbroken VT symmetry we may always change the phase 
of the initial ket vector in such a way that = In the original 
L 2 (IR) context this was a normalization convention in which both the respec- 
tive symmetric and antisymmetric components of \l/(x) = S(x) + iA(x) with 
properties S(—x) = S(x) and A(— x) = —A(x) were real. 

In the current literature, people sometimes speak about the VT symmetry 
of the system while tacitly assuming that it is not spontaneously broken, i.e., 
that the spectrum is real and non-degenerate. Under such an assumption it 
is rather straightforward to return to our specific model and to write down 
the definitions of the two new wave functions ip±{y±) with y± e (0, L) in 
terms of the components of the old wave function ^(x), 

V-(y-) = S(L-y-)-iA{L-y-), Mv+) = S(L-y + )+iA(L-y + ) . (9) 

This formula confirms that the time reversal operator T itself acts, as usual, 
as complex conjugation. 

In the new notation the general solution of differential Eq. (j3J) 

4>±(y) — ^4± sm ky + B± cos ky (10) 

must be restricted, first of all, by the external boundary conditions at y — 0. 
This yields the rule 

kA± = ±iaB±. (11) 

Its insertion in Eq. f flUj) defines the wave functions. Finally, the necessity of 

their matching in the central vertex, i.e., relation 

ktankL — ia kta.nkL + ia 
k + iat&nkL k — iatankL 
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leads to the ultimate secular equation 

k 2 - a 2 



tanfcL = 0. (13) 



k 2 + a 2 tan 2 kL 

Although this equation looks different from the secular equation as given in 
Ref. [S] (where one merely has to set d = 2L in eqs. Nr. 13 and 14), the 
set of the resulting energy roots ([3]) remains the same of course. It is worth 
noticing that from our present form of secular equation (TIB"]) the complete 
set of eigenvalues is determined via zeros of a triplet of elementary functions 
fi{k) := k 2 — a 2 , f2(k) := tan/cL and fs(k) := cotan/cL. 



3 The new model with q equilateral edges 

In the above-introduced notation the generalization of the model becomes 
straightforward. At any q > 2 we merely consider a q— plet of Schrodinger 
equations 
d 2 

- J^AVi) = E MVi) > Vi e (°> L ) > 3 = 0, 1, • • • , q ~ 1 (14) 

for which the q— plet of edges ej = ej(yj) with i/j 6 (0, L) may be visualized 
as forming a star-shaped graph with the single central vertex at yj = L. 
For the sake of simplicity, the matching in the origin will be chosen in the 
elementary Kirchhoff 's form 

1> j (L)=,/> (L), j = l,2,...,g-l, J2 9 vM L ) = °- ( 15 ) 

3=0 

In a completion of the tentative analogy, the complex rotation by angle 7r as 
used in Eq. (jBJ) will be replaced now by the complex rotation by an appro- 
priate fractional angle <fi = <p(q) = 2ir/q, yielding the prescription 

d x .iPj(0) = ict V;(0) , j = 0, 1, . . . , q - 1 , V = 2vr/g . (16) 
3.1 Wave functions 

The general solution of the differential Schrodinger system ( IT41) 

ipj(x) = Aj sin kx + Bj cos kx , j — 0, 1, . . . , q — 1 (17) 
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yields also the auxiliary expression for the derivatives, 

d x ipj(x) = k Aj cos kx — k Bj sin kx , j = 0, 1, . . . , q — 1 . (18) 



One converts the dynamical boundary conditions ( 1161) into an elementary 
connection between coefficients, 

k Aj = ia e ijtp Bj , j = 0, 1, . . . , q - 1 , V = 2ix/q. (19) 

The continuity condition for wave functions in the central vertex 

[ia e 1 ^ sin kL + k cos kL] Bj = k g , j = 0, 1, . . . , q — 1 (20) 

enables us to define all of the coefficients Bj = Bj(g, L, k) as proportional to 
the auxiliary parameter g. Their subsequent insertion in the explicit version 
?-i 

[ia e ljip cos kL — k sin kL ] Bj = (21) 

of the Kirchhoff 's law of Eq. (|15|) finally leads to a rather complicated trigono- 
metric secular equation which defines, in principle at least, all of the bound- 
state energies E = E n . As long as the underlying Hamiltonian is non- 
Hermitian, these energies may be both real and complex at q > 2. Some 
of them also need not remain expressible via any closed-form analogue of the 
special q = 2 formula (j3J). 

3.2 Secular equation 

After the elimination of BjS the secular equation for bound-state energies 
E = k 2 acquires a compactified form 

tan[kL - pj(k)] = (22) 

j=0 

where the parameter g dropped out and where we defined, implicitly, 

tan/^fc) = m eX PW) = C + iK, j = 0, 1, . . . , g — 1 , tp = 2n/q . 

Setting (3j = u + iv we may abbreviate tan u = U and tanh v — V and write 
U + iV (U + iV)(l+iUV) U{1 - V 2 ) + iV(l + U 2 ) 



tan j3j 



1-iUV 1 + U 2 V 2 ' 1 + U 2 V 2 
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As long as 

C = C(j, k) = ~ sin jtp , K = K(J, k) = y cos j<p 
k k 

we obtain the real-function correspondences 

sin u cos u sinhfcoshf 

k > = 2 , • 1,2 > K U> fc ) = — T2 , ~ 

cos z u + sinn v smh f + cos z u 

as well as the closed-form inversion formulae 

sin 2u ( C\ a 2 sinh 2 v + sin 2 u , ~ T ^ 9N 

tanjy;=— — - =- , — = , 2 — = C 2 + . 

smn 2v \ K J k z smh u + cos 2 w 

From the latter relation we may finally eliminate sinhf and convert the 

former relation into an easily solvable quadratic equation for the value of 

cos 2 u, yielding the two eligible roots as functions of j and k. 

In this manner, secular equation ( |22|) would be given a lengthy and rather 

clumsy but still explicit elementary form which we are not going to display 

here of course. Anyhow, at any pair of given parameters q and a the whole 

eigenvalue problem may be now solved, numerically, with arbitrary precision. 



4 Secular equation revisited 

Let us now demonstrate that for the purposes of symbolic-manipulation sim- 
plifications, secular Eq. (|22|) should be reconsidered in the apparently more 
complicated form of the sum 

v 1 ^ tan fcL-C(j, k) -\K(j, k) = 
1 + [CO\AO+LK-0\A:)] tanfcL ' 
In what follows we are now going to show that and how the simplification of 
this formula may be achieved by the explicit summation when proceeding, 
in a systematic inductive manner, from the smallest integers q upwards. 

4.1 The trivial single-line quantum graph: q = 2. 

At q = 2 we have tan / 9o i i(fc) = ±ia/k so that one obtains C(j, k) = and 
K({0,l},k) = ±a/k. This makes secular equation ( 1231) reducible to the 
ultimate elementary constraint ([13]) . 
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4.2 The three-pointed star graph 

When we abbreviate 7 = a/2, the q = 3 three-term secular equation 

k t&nkL — 2i7 k tan fcL + 7 v^3 + i7 tan fcL — 7 \/3 + i7 ^ 
/c + 2i7tan/cL fc - (i + ^3) 7tan kL k — (i— a/3) 7 tan /cL ~~ 

may be rewritten in its simplified, single-term form 



k 3 + ia 3 tan kL 

3 71 77T tan kL = • 25 

fc 3 -ia 3 tan 3 /cL 



The inspection of this secular equation reveals that the real and discrete part 
of the q = 3 bound-state spectrum coincides with its q = 2 predecessor, up 
to the anomalous root k = a which now disappeared. In other words, the 
real roots of the new secular Eq. fl25|) coincide now strictly with the zeros of 
the real functions tan kL and cotan kL. 

It is necessary to add that our secular Eq. f[2"5"j) also possesses complex 
roots defined by subcondition 

k 3 + ia 3 tan kL = . (26) 

We may decompose k/a = fi+iv, set A = a L and obtain the complex version 
of such a secular subequation 

3 2 . , 2 ?A tan/iA + i tanhz/A 

l — i tan /zA tanh i/A 

which is equivalent to the coupled pair of the real secular subequations 
jj 3 — 3 n v 2 + (3 /i 2 z/ — v 3 — 1/ tanh z/A) tan \i\ tanh i/A = 

and 

3 fj 2 ^ — v 3 — tanh z/A — (;U 3 — 3 /i z/ 2 ) tan /iA tanh z/A = . 

In an extensive numerical test we revealed and demonstrated the existence 
of nontrivial complex roots of these equations at the various values of A. 
For example, we localized the sample pair of roots with /x = 1.20484 and 
v = ±0.3507 at A = 1. This means that at q = 3 the non-Hermiticity of the 
Hamiltonian may probably be interpreted as "too strong". In other words, 
the Hamiltonian of the quantum version of this system cannot be Hermitized 



using just the standard techniques as reviewed in Ref. [9], at any coupling 
constant A. 

In this sense, the practical phenomenological applicability of our q = 3 
quantum graph appears restricted to the domain of non-linear optics [7] and 
to the various similar, recently popular classical-physics (or even classical- 
mechanics [19]) implementations and applications of the theory where the 
complex energies may and do find their natural physical interpretation. 

Naturally, the situation is different in quantum physics where, inside the 
physical Hilbert space H, the spectrum of any operator H representing an 
observable quantity must be real. Still, there exists a certain recently discov- 
ered [2U] space-projection trick which may prove acceptable in at least some 
phenomenological considerations and applications of our models. 

4.3 Quantum-Hilbert-space construction at q = 3 

We just demonstrated, constructively, that the discrete energy spectrum of at 
least some of our present q > 2 quantum-graph models need not necessarily 
be all real. In all of these "non-real-spectrum" cases it seems necessary to 
discard the underlying Hamiltonian H = H^^a) as leading to non-unitary 
evolution of the quantum system in question. At the same time, many of the 
formal (e.g., solvability) as well as phenomenological (e.g., scattering- related 
[21]) features of these and similar models might seem appealing enough. For 
this reason, let us now describe, briefly, one of the recently discovered and 
more or less universal remedies of the apparent complex-energy shortcoming. 

Firstly, let us remind the readers that within the rigorous quantum theo- 
ries one must often exclude all of the non-Hermitian Hamiltonians H (leading 
to the real spectrum or not) which cannot be assigned a suitable Hilbert space 
in which they may be reinterpreted (typically, via a suitable inner product 
[21 [9]) as self-adjoint. In particular, this year it has been proved [22] that in 
this manner it would be even necessary to discard the popular imaginary cu- 
bic oscillator and many other standard benchmark VT— symmetric quantum 
models with real spectra. 
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Naturally, this conclusion may seem rather surprising. Unfortunately, it is 
based on the rigorous functional analysis and, hence, mathematically valid. 
The essence of the apparent paradox has been found in the inconsistency 
between the a priori choice of the same domain T> of H before and after the 
Hermitization, V^ ciotc \H) = £>( aftcr )(#). 

This conclusion may be perceived as indication of the way out of the 
trap. Indeed, a physics-oriented and pragmatic (thought still mathematically 
rigorous) way out of such a form of crisis of the theory has been found, 
almost in parallel, in Ref. [20]. In a simplified explanation it has been merely 
admitted that T>( heioIe \H) ^ V^ {tct \H). 

The resulting flexibility of the "projection" on the meaningful vector space 
X)( after ) (iJ) of the correct physical states enables us to construct the latter 
space simply as spanned by any subset of the eigenvectors of the Hamiltonian 
H in question. In other words, the formal recipe as presented in Ref. J2U] may 
be simply read as just another version of the innovative implementation of 
the abstract principles of quantum theory where the correct Hilbert space % 
is determined dynamically (plus, in the present case, in the mere real-energy 
subspace) . 

Once we return now to our present specific q = 3 quantum-graph Hamil- 
tonian H where we choose, a priori, the usual and friendly (but, in general, 
unphysical) Hilbert-space domain t>Q**<™)(h) = 0j =o L 2 (ej) (in this space, 
H iff of course), we may now use just the slightly adapted recipe of 
Ref. [20]. Thus, in essence, we have to construct the correct vector space 
X>( after ) (77) (as well as its bra- vector dual) as spanned just by the real- 
eigenvalue eigenvectors of H (or of W, respectively). 

We omit the further details here, summarizing that due to the infinite 
number of the real eigenvalues at our disposal, the dimension of T>( after )(if ) 
will remain infinite. The immanent projector-operator nature of the whole 
construction may, indeed, be perceived as far from trivial, with details lying, 
certainly, far beyond of the scope of our present paper. Thus, we may only 
add that although, in the projected-space approach, the subsequent Hermi- 
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tization of the model remains entirely routine [9] , the physical interpretation 
or our star-shaped-graph models remains the same as in the most elementary 
q = 2 special case. For this reason, some of the apparent paradoxes (like, 
e.g., the intrinsically non-local nature of the well known q = 2 system |21j ) 
will survive the transition to q > 2 of course. 



4.4 Secular equation at q = 4 

The q = 4 version of our secular equation reads 

k tan kL — ia k tan kL + ia k tan kL — a k tan kL + a 



+ 



+ 



+ 



k + i a tan kL k — iatankL k + at&nkL k— a tan kL 

and may be again simplified, 

k 4 + a 4 tan 2 kL 

4 ~t~a it, — tan kL = . 

k 4 - a 4 tan 4 kL 



(27) 



(28) 



The real stable-bound-state spectrum remains the same as at q = 3. The 
complex roots of the auxiliary subequation k 4 + a 4 tan 2 kL = 0, i.e., of the 
two equations k 2 = ±ia 2 tankL related by the formal change of k — > —k 
may be sought just in a half-plane of complex k. The final analysis of this 
equation may again proceed in the manner outlined in preceding subsection. 




Figure 1: Numerical identification of a complex root of secular equation. 

Once we make a particular choice of the sign (say, plus), we obtain the 
simplest special case of the equation for the complex roots k = fi + i v. In 
units L = 1 we have 



/i — v +2i/!^ = +ia: tan(/i + i v) 



(29) 



12 



This is a complex equation which decays into the pair of real conditions 

p 2 — v 2 + 2 pv tan p tanh v + a 2 tanh v = , (30) 

2pv — (/i 2 — i/ 2 ) tan /i tanh z/ — a 2 tan /i = . (31) 

The search for a nontrivial root of these two equations remains numerical. 
Figure [T] samples the localization of a complex root k = \/E with components 
fi sa 1.7025 and i/ « -0.3165 at a = 1. 

Let us add that the efficiency as well as the reliability of the latter search 
for a complex root was significantly enhanced by a trick in which the two real 
equations (130]) and fl3Tl) (defining, in our illustrative picture, curves A and B 
in the p — v plane, respectively) were complemented by the third, redundant 
rule of the coincidence of the absolute values of the left- and right-hand side 
of the complex relation k 2 = ±ia 2 tan k, 



> 2 2 tan p + tanh v 



^ + ^ = a\ r 2 . (32) 

1 + tan jj tanh v 

In our picture this defines the redundant, additional curve C in the \x — v 
plane. We see that its construction helps us to identify the root in question 
via a triple intersection of the curves A, B and C. Naturally, such a recipe 
keeps the numerical errors under a very reliable control. 

4.5 Secular equation at q = 5 

At q = 5 our secular equation becomes perceivably more complicated but the 
use of computerized symbolic manipulations and appropriate trigonometric 
identities is still found to lead to the thoroughly simplified prescription 

k 5 — ict 5 tan 3 kL 

tanA;L = 0, q — o. (33) 



k 5 + ia 5 tan 5 kL 

This formula indicates that the real part of the spectrum remains unchanged 
also at q — 5. On the more important methodical level the structure of this 
formula confirms our expectation that the sum (|23|) may be represented by 
a very simple function of q. 
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4.6 Secular equation at any q 

On the basis of the above particular results yielding the explicit and ele- 
mentary closed summation formulae it is straightforward to conjecture and 
prove the validity of the following extrapolated general trigonometric secular 
equation 

k q + (ia) q t&n q - 2 kL , rn 

Whenever we restrict our attention just to the real roots k which correspond 
to the stable bound-state solutions, we reveal that the numerator in the 
fractional part of the left-hand-side secular determinant (134")) plays now the 
role of the source of the q > 2 analogues of the single anomalous q = 2 root 
k = a. 

This numerator cannot vanish at real k and odd q = 2m + 1, m G Z + and 
it cannot vanish at q — 4m, m G Z + , either. The remaining values of the 
integer q = 4m — 2, m G Z + appear exceptional. Their choice leads to the 
emergence of the additional real zeros and so it deserves a separate attention. 

4.7 Secular equation at the exceptional q = 4m — 2 

At q = 2 our quantum-graph spectrum was all real (cf. sec. 14.11 above) 
but, as we saw, the situation became anomalous at any q > 3. Nevertheless, 
what is new at the exceptional integers q = 4m — 2, m = 2, 3, ... is that the 
secular-equation factor 

k Am-2 _ a Am-2 ^m-A fc£ _ Q (35) 

becomes nontrivial and, moreover, that it produces, obviously, certain po- 
tentially real additional bound-state eigenvalues E = k 2 . 

One of the unfortunate consequences of the latter observation is that some 
of the energies of the stable quantum-star bound states cease to be obtainable 
in closed form. Their determination must be performed by suitable brute- 
force numerical methods. Moreover, the practical search for the roots of 
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transcendental Eq. (1351) or of its slightly more friendly and graphically better 
tractable version 



^anomalous Oi 



tan i-i/( 2m -i) kanomalousL = o , m = 1, 2, . . . (36) 



becomes technically complicated. A priori, without any extensive calcula- 
tions we can immediately be sure that at the sufficiently small as, all of the 
generic and m— independent real roots k — (n + l/2)n ^> 1 with n — 0, 1, . . . 
become accompanied by the neighboring real pairs of eigenvalues which are 
produced by Eq. f l3"6"|) . 

With the growth of a this picture will first lose its validity between k = 
and k = tt/2. At the sufficiently small as one always finds there the two 
smallest real roots inside the interval. With the growth of parameter a these 
two roots move towards each other at a speed which depends on m. In a 
numerical experiment performed at m = 2 we found that there exists the 
critical value of a = a cr iu C ai ~ 0.7863 at which these two lowest anomalous 
energy- level twins merge at k « 0.748 and, subsequently, complexify. 

This observation may only be read as a reliable numerical proof that at 
the sufficiently large values of the strength of the non-Hermiticity a > a cr i tica i 
(with ot cr iti ca i ~ 0.7863 at m = 2), the spectrum of the whole system certainly 
contains non-real eigenvalues. 

In the interval of a < 0.7863 we may only conclude that there exists a set 
of certain new and strictly real "anomalous" eigenvalues which may only be 
generated numerically (i.e., say, via our secular sub-equation (I3"6"|) ). This ex- 
tremely interesting infinite family of the new quantum states is, in principle, 
observable. Its energy levels (which form, incidentally, almost degenerate 
pairs at higher excitations) may be interpreted as the appropriate quantum- 
graph q = 6, q = 10 (etc) analogues of their single-state q = 2 predecessor 
E (a) of Eq. ©. 
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5 Summary 



In our preceding papers [23j El] on VT— symmetric quantum graphs we al- 
ways restricted our attention to their mere discrete approximants. In our 
present paper, we abandoned this approach as not sufficiently efficient. An 
alternative way of circumventing the technical obstacles has been found here 
in an ad hoc restriction of the class of the admissible graphs to their star- 
shaped subset G^ q \ Due to this restriction, we were able to replace the 
universal though less powerful discretization approach by the much more 
elementary method of matching of wave functions at the central vertex. 

In technical sense, our present results may be perceived as a return to 
optimism. The main source of the simplification of our constructive con- 
siderations may be identified with the inherent symmetry of the complex 
Robin boundary conditions. This symmetry found its fructification in the 
emergence of powerful trigonometric identities. These identities led to the 
enormous simplification of the related secular equations at any integer q > 2. 
In this sense one could find here certain parallels with the role of trigonomet- 
ric identities, say, during the early stages of development of Calogero models 
[25] and/or of some of their less influential analogues [26J. 

The transition to nontrivial topology of the graph-related phase space 
(or of the space of coordinates) manifested itself in two ways. Firstly, a part 
of the physical sector where the bound-state energies remained strictly real 
appeared independent of the number of rays q, i.e., mathematically stable. 
Secondly, strictly this part of the spectrum also remained defined by closed 
formulae at q > 2. In contrast, the rest of the spectrum (and, in particular, 
the whole sector of "resonances" where the energies are complex) appeared 
changing with the changes of q. 

Due to the elementary form of the secular equations at any q, the "friendly" 
closed-form real energies coincided with their elementary square-well q = 2 
predecessors. In contrast, it appeared rather difficult to localize the pre- 
cise position of all of the complex bound-state energies in complex plane. A 
sophisticated numerical approach appeared necessary for the purpose. 
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In the context of physics the potential phenomenological applicability of 
our present family of toy-model quantum graphs may be perceived as guided 
by the parallels with the q = 2 special-case square-well which represents one 
of the simplest available VT— symmetric models. This parallelism may be 
expected to include, e.g., a potential relation between the present q = 2 non- 
constant level Eq = Eq{q) of Eq. (|3]) and the similar anomalous levels which 
are known to emerge in supersymmetric models [27] . In the future, other pos- 
sible parallels might also appear reflecting, say, the preservation of a certain 
complex- rotational symmetry of our present wave functions (cf. Eq. ([16]) ) or 
the related graph-inspired permutation-transformation generalization of the 
concept of the parity, etc. 

In the context of mathematics, one of the most unexpected byproducts of 
the transition to q > 2 occurred at the subsequence of models with q = Am— 2. 
In contrast to the presence of a single anomalous real energy level with k — a 
which existed at m = 1, it has been found that infinitely many anomalous 
real energy levels seem to exist at any larger m > 2. This phenomenon is a 
truly puzzling new structural feature of the spectrum of a phenomenological 
model. Its deeper theoretical explanation (say, via its possible relation to the 
complex-rotational symmetries of wave functions) remains an open question 
at present. 
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